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Abstract 

The collapse of a rotating 3-axis ellipsoid is approximated by a system of ordi- 
nary differential equations. Violent relaxation, mass and angular momentum 
losses are taken into account phenomenologically. The formation of the equi- 
librium configuration and different types of instability are investigated. 
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According to modern cosmological ideas the most part of matter in the Universe 
belongs to so called Cold Dark Matter, consisting of non-relativistic particles. 

The study of the formation of dark matter objects in the Universe is based on 
N-body simulations, which are very time consuming. In this situation a simplified 
approach may become useful. 

Here we derive equations for the dynamical behaviour of a compressible rotating 
3-axis ellipsoid in which a motion along axes takes place in the common gravita- 
tional field of a uniform ellipsoid and under the action of the isotropic pressure. 
Collapse in the dark matter are characterized by non-collisional relaxation, based 
on the idea of a "violent relaxation" of Lynden-Bell PJ. So the collapse of the el- 
lipsoid is approximated by a system of ordinary differential equations, where the 
relaxation and losses of energy, mass and angular momentum are taken into account 
phenomenologically. The system is solved numerically. 

Bisnovatyi-Kogan 1 considered the case of spheroid (a = b ^ c), where there 
are analytical formulaes for the gravitational potential and forces. 

Let us consider a uniform 3-axis ellipsoid, consisted of non-collisional non- 
relativistic particles, with semi-axes a =/= b =/= c and rotating uniformly with an 
angular velocity f2 around the axis z. 

The mass m and total angular momentum M of a uniform ellipsoid are connected 
with uniform density, angular velocity and semi-axes as m = 4^ p abc , M = 
Y ft(a 2 + b 2 ) . Assume a linear dependence of the velocity on the coordinates: v x = 
ax I a , v y = by/b , v z = cz/c . The gravitational energy of the uniform ellipsoid is 
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defined as: 

oc 

_ 3Gto 2 r du 

9 ~ ItT J ^(a 2 +u)(b 2 +u)(c 2 + U ) 

and is expressed in elliptical integrals. 

Consider a compressible ellipsoid with a constant mass and angular momentum, 
a total thermal energy of non-relativistic dark matter particles E t h ~ V~ 2 ' 3 ~ 
(a&c) -2 / 3 and the relation between pressure P and thermal energy E t h as Eth = %y- 
In absence of any dissipation the ellipsoid is a conservative system. 

To derive equations of motion let write for it a Lagrange function 

L = Ukin - U p ot , U p ot = U g + Eth + Urot , Ukin = yjj + b 2 + (?) , 
c E th ^ n (a in b in Ci n ) 2 / 3 £ rr 1 /" 2 5 M 2 
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By variation of the Lagrange function we obtain Lagrange equations of motion. It is 
easy to check that equilibrium solution of these equations is the Maclaurin spheroid, 
and in addition the 3-axis Jacobi ellipsoid appears at large angular momentum. 

In reality there is the "violent relaxation" in the collisionless system 0. There- 
fore there is a drag force, which is described phenomenologically by adding of the 
terms — 7^7 , — 7^7 , — 7^7 in the right-hand parts of equations of motions. 

Here we have scaled the relaxation time r re i by the Jeans characteristic time 
with a constant value of a re i : r re i = a Te \Tj — 2w a re i . 

The process of relaxation is accompanied also by energy, mass and angular mo- 
mentum losses from the system. These losses may be described phenomenologically 
by characteristic times r e ; , r m ; , tmi ■ 

To obtain a numerical solution of equations we write them using non-dimensional 
variables. Taking into account mass and angular momentum losses, the dynamics 
of the system is described by the following non-dimensional system of equations 
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m dt 2 J (a 2 + u)y/(a 2 + u)(b 2 + u)(c 2 + u) 3m a (abc) 2 / 3 
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00 

■■ b dm 3m f du Idle 



111 



dt 2 J (&2 + u )yJ{a 2 + u){b 2 + u)(c 2 + u) 3m b (abc) 2 / 3 

25M 2 b b 
m 2 (a 2 + 6 2 ) 2 T re i ' 
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„ c dm 3m f du in 1 s 
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dt 2 J ( c 2 + u ) V(a 2 + u)(b 2 + n)(c 2 + u) 3m c (abc) 2 / 3 T rel 
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This system is solved numerically for several initial parameters. 

The case of spheroid (a = b ^ c) is considered by Bisnovatyi-Kogan £Q. Our 3-d 
equations give the same result for spheroidal initial conditions. 

In the case of 3-axis ellipsoid there are new qualitative effects, because there 
is additional degree of freedom compared to the case of spheroid. At large initial 
angular momentum Mi n we observe the instability of the Maclaurin spheroid for 
the transformation into the Jacobi ellipsoid. So in this case the system reaches 
equilibrium 3-axis ellipsoid. See Fig.l. 




Figure 1: The instability at large momentum. 

At small initial angular momentum Mi n the system reaches the oblate spheroid. 
But we have an instability, disappeared during the relaxation. It is the instability 
of system with purely radial trajectories. 

The work is in progress. We are going to investigate the different aspects of 
instability in more detail. 
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